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Abstract 

This paper presents some algorithmic techniques to compute exphcitly the noetherian oper- 
ators associated to a class of ideals and modules over a polynomial ring. The procedures we 
include in this work can be easily encoded in computer algebra packages such as CoCoA [5] . 

1 Introduction 

The Ehrenpreis-Palamodov Fundamental Principle, [7j and [19], states the following: 

Theorem 1.1. Let pi{D), . . . ,pr{D) be linear constant coefficients partial differential operators 
in n variables. Then there are algebraic varieties Vi,...,Vf in C" and differential operators 
di, . . . ,dt with polynomial coefficients, such that every function f G C°°(M") satisfying 

p,{D)f = ...=pr{D)f = 

can be represented as 

f{x) = j2[ 9,(e"-)^^,(z), (1) 

for suitable Radon measures duj . 
The collection 

V = {{Vi,diy,{V2,d2);...;{Vt,dt)} 

is said to be a multiplicity variety and Theorem 11.11 is equivalent to the following strengthening 
of the classic Nullstellensatz: 

Theorem 1.2. Let L be an ideal of R. There exists a multiplicity variety V such that a poly- 
nomial f belongs to I if and only if djf\y. = for every j = 1, . . . , t. 
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The operators di,...,dt are called, in Palamodov's terminology, noetherian operators be- 
cause their construction relies essentially on a theorem of M. Noether on a membership criterium 
for polynomial submodules (see e.g. [19] pp.161, 162). The nature of the original proof of the 
Fundamental Principle is essentially existential and therefore the question of the explicit con- 
struction of such operators is of great interest whenever we consider a concrete application of 
the Fundamental Principle. Note that if / is the ideal generated by the polynomials pi, . . . ,pr 
and if 

I = n---n Qt 

is its primary decomposition, then the varieties Vj which appear in theorem 1 1.1 1 are simply given 
by the algebraic sets V{Qj). The information on the multiplicity of each of them is left to the 
operators dj. 

In this paper we build on some recent results in the construction of noetherian operators 
[TI] ITSl ITTl ITSl [20] and provide some new algorithms which allow the automatic construction 
of these operators at least in some rather large class of cases. We include several experiments 
using algorithms implemented on CoCoA 

In section 2 we quickly review the fundamental tools from computational algebra (mostly 
the theory of Grobner Bases). The core of the paper is section 3 where we deal with case of 
zerodimensional ideals and where we present several explicit algorithms. A final section deals 
with the case of ideals of positive dimension. 

Executable versions of the algorithms discussed in this paper have been explicitly written 
for CoCqA. and are freely available at 

\protect\vrule widthOpt\protect\href {http: //www. tlcl85 . com/ coala}{http : //www. tlcl85 . com/ coala^ 

The authors are grateful to F. Colombo, S. Hosten, and B. Sturmfels for their many useful 
suggestions and comments. The first author is grateful to George Mason University and to the 
Eduard Cech Center for the financial support. 

2 Computational Algebra Tools 

Throughout this paper, we will work in the ring R = C[xi, . . . , Xn] of polynomials in n variables 
with complex coefficients; we will think of R as the ring of symbols for the differential operators 
we are studying. Even though we consider differential operators with constant coefficients, the 
Fundamental Principle shows that noetherian operators have, in general, polynomial coefficients; 
we will use the symbol An to denote the Weyl Algebra C[xi, . . . , Xn, dxi, . . . , dxn] of such oper- 
ators. Here, and throughout the paper, the symbol dx will be used as a shortcut for 
Using the notation introduced in [12], we will denote the monoid of power products in R by T" 
and the module monoid of power products in R^ by 

T"(ei,...,e,) = {iei | t G T", i = 1, . . . , s} 

where Cj is the i-th element of the canonical basis of R^. All the definitions to follow are given 
for ideals but can be extended in straightforward fashion to the case of modules |12J. A term 
ordering a on T" is a total ordering on power products with the following two properties: 

I) if h t2 and t £ T" then t-ti>„t- 12; 

II) if t e T" and s G T" then s-t>^t. 
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The leading term ideal associated to / with respect to a is the ideal generated by all the leading 
terms of elements of /, and will be indicated by 

LT.(/) = ({LT,(/)|/G/}). 

More in general, the leading term ideal associated to a subset G oi R will be written as LTo-(G) = 
({LT„(/)|/ e G}). Note that LT„{G) = LT„{I) if and only if the set G is a Grobner Basis for 
the ideal /, this being the main characterization of a Grobner Basis. 

The algorithm which associates to an ideal I oi R its Grobner Basis Ga{I) is the core algorithm 
of the theory of Grobner Bases and can be found for example in [12j . theorem 2.5.5. Another 
key tool in computational algebra is the division algorithm (see again |12) . theorem 1.6.4) which 
can be performed to generate the remainder of a polynomial with respect to a set of generators 
of /. Note that the remainder of a polynomial depends on the set of generators chosen for / 
(in fact, it even depends on their order). The fundamental property of Grobner Bases is that 
such a remainder is zero if and only if the polynomial belongs to the ideal. For this reason the 
remainder calculated with respect to a Grobner Basis is called the normal form of a polynomial. 
Given a polynomial f € I and a term ordering a, we will denote by NFg- (/) the normal form 
of / with respect to the cr-Grobner Basis of / (the same notation is used for modules). An 
equivalent way to compute a remainder is using rewrite rules (see |12j section 2.2). Given a 
polynomial g €z R, we say that a polynomial /i rewrites to /2 with respect to the rewrite rule 
— ^ (and this is indicated by /i — ^ /2) if there exists a monomial m in i? such that /2 = /i — rng 
and LT f^{mg) is not in the support of f2- This is also called a one-step reduction. We can rewrite 
a polynomial using a set of elements Q = {gi, . . . ,gs} by performing a one-step reduction with 

each of the gi^s, in that order. We will denote by — > the transitive closure of the relations 
. . . , This relation is called rewrite relation or rewrite rule. By applying a sequence of 
one-step reductions to a polynomial / using the elements in G we then obtain a remainder of / 
with respect to {gi, . . . ,gs}- In particular if is a Grobner Basis we have that / rewrites to its 

normal form, i.e. /-^NFo- (/). 

We now introduce some definitions about elimination theory and term orderings (see |12j . 
section 3.4, for details on this topic) which will be necessary in the last section. 

Definition 2.1. Let R = C[x,t] where x = (xi, . . . ,Xn-d), t = (ti, . . . ,td). A term ordering a 
on T" is called an elimination ordering with respect to x if every element f £ R whose leding 
term is contained in C[t] is such that / G C[t]. In other words, 

V/ G R, LTM) G m ^ / G C[t]. 

The reason why such a term ordering is called an elimination ordering is that it allows to 
eliminate the variables x from an ideal, i.e. it allows to compute ID C[t]. To do this, it suffices 
to compute a Grobner Basis with respect to any elimination ordering as in definition 12.11 and 
then keep only the elements that do not contain any monomials in x. Such elements actually 
form a Grobner Basis for the ideal In C[t]. It can be easily checked that Lex, the lexicographic 
term ordering on T", is an elimination ordering with respect to any "initial" subset of variables, 
i.e. with respect to any subset of the type {xi, . . . , Xk} in C[xi, . . . , with k < n. A class of 
term orderings that satisfy the elimination property and that we are going to use for our goal 
of computing the noetherian operators in C(t)[x] are the so called product orderings. 

Definition 2.2. Let R = C[x,t] as before and let ax and at be two term orderings on the set of 
terms Tx = {x"' \ a G N"~'^} and Tt = {t'' \ b G N'^} respectively. The product ordering ax • at is 
defined by 

xH'' >ax-at xH'^ ^ x" or (x" = x" andt^ t"^). 
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It is immediate to show that the product ordering defined above is an ehmination ordering 
with respect to x, no matter what the choice of ax and at is. Ehmination orderings are usually 
slow when it comes to Grobner Basis computations, in particularly Lex is known to be one of 
the slowest. Product orderings are then introduced to perform better. One can in fact define a 
"fast" term ordering (such as DegRevLex) on each of the two subsets of variables, and then take 
the product. The following lemma will be useful later in the paper. 

Lemma 2.3. Let R = C[x,t] he a polynomial ring equipped with a product ordering a of the 
type ax ■ at as in definition \ 2. S[ Let I he an ideal of R and let Q = {gi, . . . , gg) he a cr -Grobner 
Basis for I. Consider the extended ideal IRd in Rd = C{t)[x\ endowed with the term ordering 
ax ■ Then Q forms a Grobner Basis for IR^ with respect to ax ■ 

Proof. Denote with x"'H'^^ the leading term of gi, where a, G and Cj G N'^, i = 1, . . . , s. 

From the fact that we chose a product ordering a, it follows that once we view gi as an element 
of IRd, its leading term is x"'. In other words, UTa^idi) = a;"* in Rd- Consider a polynomial / 
in IRd- The set Q still forms a set of generators for the extended ideal, so / can be written as 
an i?(i-linear combination of the giS. Moreover, supposing / monic, we can write / as 

f = x"" + ^p^{t)x^ , where 5 G N'^-"' and x^Mh. 

b 

Consider the product D{t) of all the denominators of the coefficients Pb{t) in /. Then D{t)f 
is a polynomial in R and it is still a combination of the elements of so D{t)f G /. Because 
of the fact that cr is a product order, the leading term of D{t)f is simply the leading term of 
/ multiplied by some power of t, i.e. lJT„{D{t)f) = xH" for some c G N^. Hence, G bemg a 
Grobner Basis for /, x"'t'^ is a multiple of one of the leading terms of its elements, say x^-^t^^ 
modulo a change on the order in Q. This means that there exist a G and 7 G N*^ such that 

x^t^ = x'^t^x'^H"^ 

which means that is a multiple of x"'^ , and this concludes the proof. □ 

3 The Zerodimensional case 

In this section, / is a primary zerodimensional ideal, i.e. the algebraic set V(/) is a finite union 
of points in C". Since a zerodimensional primary ideal is associated to a single point of the 
variety V(/) we can always assume, with a change of coordinates, that V(/) = {(0, . . . ,0)}, or 
equivalently that vT = (xi, . . . , 

3.1 Closed Differential Conditions 

A first complete description of the differential condition characterizing a zerodimensional primary 
ideal centered in zero has been done in [15]: we briefly recall the main notations and definitions 
of that paper. We will denote with D{i\, . . . ,in ) : ii — > i? the differential operator defined by: 

D{ii,...,in) = . , V . dx\' ■ ■ ■ dx'^ , G N , forahj = 1,... ,n, 

III - ■ -Inl 

or, alternatively, if i = x^^ • • • x^" G T", we will use the symbol D{t) as D{ii, . . . , in). Moreover, 
we write V = {D{t)\t G T"} and denote by Spanc(I') the C- vector space generated by V. We 
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now introduce some morphisms on T) that act as "derivative" and "integral": 

[0 otiierwise 

Px^ {D{ii, . . .,in)) = D{ii, . . . ,ij + 1, . . . (3) 

Such operators extend trivially on Spanc(I') by linearity, and one can easily define at and pt 
for any t G T" by composition. 

Definition 3.1. A subspace L of Spanc{T>) is said to be closed if 

O'xj {L) ^ L, for all j = 1, . . . , n. 

Definition 3.2. Let I be a primary ideal in R such that \fl = {xi, . . . ,Xn)- We define the 
subspace of differential operators associated to I as 

A{I) := {LG Spanc(P)|L(/)(0,...,0) =0/ora///Gl}. 

Similarly, we associate to each subset V C Spanc(I') an ideal 

I{V) :={feR\ L(/)(0, . . . , 0) = 0/or all L e V} 

Theorem 3.3. Let m be the maximal ideal {x\, . . . , Xn) of R. There is a bijective correspondence 
between m-primary ideals of R and closed subspaces of SpanciV) 

A 

{m — primary ideals in R} T^{closed subspaces of Span^{T>)} 

I 

so that I = ZA{V) and V = AI{I) for every I and V. Moreover, for a zerodimensional 
m-primary ideal of R whose multiplicity is p, we have that dimc(A(/)) = p. 

Theorem 13.31 shows that the noetherian operators associated to a zerodimensional primary 
ideal form a closed subspace of Spanc(2?). In addition, when considering a zerodimensional 
primary ideal, since the dimension of A(/) is finite, we can view a basis of A(/) as a set of 
noetherian operators which, in this particular case, happen to be operators with constant co- 
efficients. Moreover, such a vector space has the nice property of being closed, fact that has 
been used by the authors of to construct a procedure that, given /, computes A(/). The 
algorithm is described below. 

Algorithm 3.4. Let I be a zerodimensional primary ideal of R such that V{I) = {(0, . . . , 0)}and 
let p = dimc{R/I) be its multiplicity. The following procedure computes the noetherian opera- 
tors associated to I: 

Input; Q = {gi, . . . , gt} a Grobner Basis for L 
Output.- A{I) = {Lo,...,L^_i} 
Initialization; i = I, Lq = 1 =Idspa.nc{V) 

If p > 1, construct a linear operator Li = X]j=i ^j^^i ^''^ opportune choice 
of the Cj 's such that Li(/)(0, ... ,0) = is satisfied for each generator f of L 
Put i = 2. 
While i < p do 
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define Lj+i as a linear combination of pj(,(Lo)) • • • ^ Pjiil^i) such that 

- {Lq, . . . , -Li+i) is closed and 

- Lj_|_i(/)(0) = for each generator f of I 

Corollary 3.5. Let L be an operator of A{I), where I is as in algorithm \3.4\ and is its 
multiplicity. Then deg(L) < ^ as an element of An- 

Proof. The construction of A(/) starts with Lq = 1 and at each step the degree of Lj+i increases 
of at most 1, so that the last element L^_i has degree at most /i. □ 

Remark 3.6. Algorithm 13.41 consists basically in the solution of a system of linear equations in 
the coefficients Cj of the linear combinations Li+i = coPJoLq + • • • + Cipj-Li. Since the system 
can have more than one solution, one may simply pick the one with minimal norm. An imple- 
mentation for a simplified version of 13. 41 has been coded for CoCoA and is available through the 
CoAL^ webpage [1] . 

Example 3.7. The following example is taken from [7] (p. 37, ex. 4). Here we show how to study 
it using algorithm 13. 4[ Let us consider the primary ideal at the origin / = (y^, — y) C C[x, y] 
whose multiplicity is 4. We start with Lq = 1 and an obvious choice for a linear operator is 
Li = dx. This has also a geometric interpretation: the origin is the intersection of the two 
curves given by the generators y"^ (the x-axis twice) and — y parabola with vertex at the 
origin). Such two curves not only intersect at the origin but they are also tangent along the 
direction of the x-axis, therefore Li = dx must be a noetherian operator. The higher degree 
operators describe a higher contact of the line and the parabola at zero. We can try to find the 
next one as a combination L2 = adx + bdxy. However, this operator L2 does not respect the 
closure condition since (Tx{L2) = a + bdy which is not in the subspace {Lq,Li) = {l,dx). A 
different choice for the morphisms p^j , instead, gives L2 = apy{l) + bpx{dx) = ady + bdx"^ which 
respects closure and annihilates the generators of I at zero with a = 1 and b = i^- Again, this 
operator could have been foreseen in advance since it is the global annihilator of — y and it 
annihilates at the origin. As a last operator, one can choose L3 = Px{L2) = dxy + ^dx^. Of 
course, the choice Py{L2) = \dy'^ + ^dx^y would have been possible as far as the annihilation 
of I is concerned, but it would have violated closure since ax{L2) = dxy is not a combination of 
the previous operators. The iteration ends here since we have found 4 differential operators. 

3.2 Forward reduction 

We are now going to present an alternative procedure to compute the noetherian operators 
associated to / that makes no use of linear algebra and utilizes the power of Grobner Bases. 

Algorithm 3.8 (Computation of noetherian operators for zero dimensional ideals). Let L be 
a zerodimensional primary ideal of R such that V{I) = {(0, . . . ,0)}. The following procedure 
computes the noetherian operators associated to I: 

Input; Q = {gi, . . . , gt} a Grobner Basis for L. 
Output; A(/) ={Li,...,L^}. 
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• Compute fj,{I) = dimc(i?//). 

• Write the Taylor expansion at the origin of a polynomial h € R 
up to the degree fj, — 1 with coefficients G C; 

T^^\h{x\ , . . . , Xji) — ^ ^ CaX^ . . . Xy^ 

aeN", \a\<n 

• Write the Normal Form ofT^^^ih with respect to Q as 

NF^ T^-Mxi, ...,Xn) = Yl df^xf' ...xt (4) 

and find scalars a^a G C such that dp = a^aCa- 

• For each (3 such that dp / 0, return the operator 

Lf3 = ^ apa — 7-^ -.dx'l^ ■ ■ ■ ^a;"" = ^ a^„Z?(ai . . . a„). 

ai. ^ 

Proof. Let h{xi, . . . , Xn) = Yl'\a\lS Cq,x"^ . . . x"" be the Taylor expansion centered at the origin 
of a polynomial h G R and let y be the Grobner Basis of I. From the theory of Grobner Bases 
we know that the normal form with respect to ^ of /i is zero if and only if h £ I, so the condition 
NFo- h = is the one that we want to characterize. It suffices to write 

degih) deg(h) 

NF. ( j2 ca^r ...<")= x; . . . = (5) 

|a|=0 \(3\=0 

and deduce from the annihilation of each coefficient dp in ^ a differential condition on h. This 
completely characterizes the membership of a polynomial h to /. The only thing to observe is 
that we do not need to work with terms up to deg{h) for the Taylor expansion. In fact, the 
number of differential conditions we need is precisely n, and so from corollary 13.51 it follows that 
the derivatives to be considered are, in the worst case, the ones of order ^ — 1 (see also [16j). 
Those differential conditions arise by using coefficients up to \a\ = fi — 1. Therefore the 
Taylor expansion can be truncated at ^ — 1. □ 

Remark 3.9. It is crucial to observe that we do not need to characterize the membership of a 
polynomial h of undetermined degree deg{h) since we have the bound ;U — 1 on its degree. Thus 
algorithm 13.81 is a procedure that is implementable on any computer algebra software package. 
Moreover, the computation of the normal form can be done degree by degree, so that we 
can stop the reduction process whenever the normal form of a particular degree is zero. This 
actually speeds up the computations in most cases we studied (up to date CPU times for several 
example are available on ^). 

Example 3.10. Consider again Example 13.71 to show the substantial difference between pro- 
cedures 13.41 and 13. 8i Since /u(I) = 4 we start by writing the truncated Taylor expansion of a 
polynomial h € C[a;, y]: 

T3h{x, y) = Coo + ciQX + coiy + C2qx^ + cnxy + co2y^ + csoa;^ + C2ix^y + c^xy"^ + cosy^ 
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and perform the normal form computation using ^ y and ^ as rewrite rules. Grouping 
like terms we can write the remainder of T^h as a linear combination of the generators l,x,y,xy 
of R/I as follows: 

NFct (Ts/i) = [coo] + [cio]x + [coi + C2o]y + [cu + C3o]xy (6) 

We call these four terms a Macaulay basis for the ideal I, although this name is also used by 
some authors for a generalization of a Grobner Basis. Note that the terms y"^ , x'^y, xy'^ and 
disappeared since they all rewrote to zero. The computation ends by expressing the coefficients 
written into square brackets in ([6]) as operators according to their meaning as Taylor coefficients. 
Namely [cqo] 1, [cio] dx, [cqi + C20] dy + \dx'^ , [cn + C30] — > dxy + \dx^ . This gives the 
same result obtained in the example 13.71 as expected. This is not surprising since theorem 13.31 
states that the correspondence I ^ ^{I) is one-to-one. 

Algorithm 13.81 does not take directly into account the closure of the space of noetherian 
operators, as algorithm 13.41 did. The fact that A(/) is closed is a general fact which follows from 
a Leibniz formula for the morphisms a^j and the fact that / is an ideal (see |14] . prop. 2.4). 
This is true not only for zerodimensional ideals but also in positive dimension, as we will see in 
section 9. We want to show that the closure of A(/) is also a direct consequence of algorithm 
13.81 and of the following property of Macaulay bases. 

Lemma 3.11. Let I G R be an ideal and letAi be the Macaulay basis of R/I, i.e. the generators 
of R/I as a C-vector space. Let Sx^ : T" — > T" be the "derivative" morphism 

Q ■ { T*^^ . . . T*^^ ^ 

Then Ai is Sj-closed for each j. 

Proof. It is known that the Macaulay basis for R/I can be computed through a Grobner Basis 
Q of /. In fact it is (see [12], theorem 1.5.7): 

M = T"\LT,(e?) 

where a is any term ordering on T". Since ^ is a Grobner Basis for /, the leading term ideal 
LTo-(/) coincides with LTct(^). Let t 7^ be a term of Ai. Suppose that there exists an 
index j such that 7^ Sj{t) ^ A4. Then Sj{t) £ LTct(^). The latter being an ideal, we have 
t = Xj ■ Sj{t) £ LTo-(^/), which is a contradiction. Note that if Sj{t) ^ M for all j, this simply 
says that t = which is again a contradiction. □ 

The morphism Sj introduced in the above lemma is the analogue of defined in section 
13. H and we will show in the next proposition that the Sj-closure of Ai is equivalent to the 
cTxj. -closure of the space of noetherian operators associated to /. 

Proposition 3.12. Let I be a zerodimensional primary ideal of R such that V(/) = {(0, . . . ,0)} 
and let O = {Lf^} be the set of operators computed with alaorithm \3.8i Then SpanciiLp}) is a 
closed subspace of Spanc (V) . 

Proof. Let Lj^ £ O, and let dp be the corresponding coefficient of the normal form NF^- (/i) as 
computed with the algorithm. Let be the term whose coefficient is dp. It is 

clear that such a term is part of the Macaulay basis of R/I since it appears in the expression of 



xl^-.-x^ '■■■xl ifi,>0 
otherwise 
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NF(j (h), which is a representation of the class of h in the quotient R/I. Denote by Fjs the set 
of operators of O such that the corresponding term in the expression of NFq- (h) divides x^: 

Fp = {L^ £ O such that x'^ I x^} 

and for each G Ffs consider = x^~'^. Since each Lj has been computed from the Taylor 
expansion of using a division algorithm that uses a Grobner Basis G of /, we have that (see |12] . 
prop. 2.2.2) if h' is such that 

x'^ = NFo- {h') and supp(/i') C supp(/i) 

then 

= t^x^ = NF^ (^)NF^ {h') = NF^ {t^h') 

i.e. the term in x^ is obtained rewriting a multiple of that part of the polynomial h which 
rewrites to x'^. By looking at the expression of is then obvious that 

since is written as a combination of Taylor coefficients corresponding to the terms of t^h' . 
It now suffices to prove that such t^'s are enough to conclude that O is closed. This is a 
consequence of the previous lemma, since all the in are associated to those terms x'^ 
of the Macaulay basis A4 that divide x^, hence from the Sj-closure of Ai we deduce that 
{x-y =st^{x^)} = {sj{x^),j = l,...,n}. □ 

3.3 Backward reduction 

We could think of performing the reduction step of the algorithm for the computation of noethe- 
rian operators for zerodimensional ideals "backwards". Instead of writing the full Taylor expan- 
sion and then using the Grobner Basis of / to rewrite it, we start from the residual monomials, 
which are easily calculated for example with CoCoA. We then "pull back" each monomial using 
the generators of / as " anti-rewrite rules". Let us explain what we mean by this. In general, 
when using a polynomial / to rewrite another polynomial g, we use its leading monomial LT(/) 
to divide the polynomials g and then we substitute each LT(/) in g with the tail of /, LT(/) — /. 
For instance, we rewrite 5 = x^ to xy using / = x^ — y, by replacing x^ in x'^ with the tail 
x^ — (x^ — v) = V- This operation, when performed using the elements of a Grobner Basis for /, 
does not alter the class of g in R/I and leads to the normal form NFo- {g). What we mean by 
" anti-rewriting" is, roughly speaking, to use the smallest monomial of /, in{f), and replace it 
with the head of the polynomial, in{f) — f. This way, from in{f) we "climb up" to find all the 
other monomials that are equivalent to in{f) modulo (/). Here is a more precise definition. 

Definition 3.13. Let f be a polynomials of R, let g he a monomial and let m = in{f) be the 
smallest term of f with respect to a given term ordering on T". We say that g rewrites backwards 
to g' in one step, using f , if m divides g and 

g = — (m - /). 

m 

Example 3.14. With this terminology, g = xy rewrites backwards to x^ using x^ — y, which 
is exactly the opposite of the standard rewrite process that leads from x^ to xy. If we use 
f = x^ + xy — 2y instead, g = xy rewrites to g' = ^x^ + |x^y. Finally, g could not be rewritten 
backwards suing x^ — since does not divide g. Notice that in general if we perform a 
one-step backward reduction and then a one-step reduction in the usual way, we obtain back g. 
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We can now apply an iteration of this procedure of rewriting backwards a monomial using a 
Grobner Basis for /. We start from a residual monomial and we rewrite it backwards using one 
generator. Then we rewrite backwards each monomial obtained after this step, if possible, using 
any element of the Grobner Basis. Technically this procedure never ends, as we can imagine, to 
obtain a new polynomial of higher degree at each step, as for example with g = x and f = — x. 
However, for the purpose of computing noetherian operators, we know from section [3] that, as 
polynomials in C[dxi, . . . ,dxn], they have degree at most fi — 1. Therefore we can stop the 
iteration once we have reached a polynomial of such a degree. Let us illustrate this idea with 
an example before we present the algorithm in general. 

Example 3.15. Consider the ideal J = (x'^ — z, y'^ — z, z^) in C[x, y, z\. It represents the origin in 
with multiplicity eight. Its generators are a DegLex Grobner Basis. The residual monomials 
for R/J are 

{l,x,y,xy,z,xz,yz,xyz}. 

First, let us reconstruct the noetherian operators associated to xyz. By rewriting it using 
— z we obtain the new monomial x^y. This cannot be rewritten further. However, the 
term xy^ is another monomial that is "attracted" by xyz via the other generator y"^ — z oi J. 
Summing up the residual monomial and all the results of the backward reduction we then obtain 
g' = x^y + xy^ + xyz whose dual D{g') = ^dx^dy + ^dxdy^ + dxdydz is actually the noetherian 
operator of J relative to xyz. 

The choice of the residual monomial xyz in Example l3.15l is not random. Indeed it is maximal 
among all the residual monomials with respect to the derivative morphisms ([7]). 

Definition 3.16. Let m be a residual monomial of R/I. We say that m is a corner monomial 
if it is maximal with respect to the monoid structure o/T", i.e. if 

Xi ■ m £ LT(/), for all i = 1 . . . n. 

If we represent R/I as a subset of N", the corner monomials are exactly in corner position. 
Proposition 13. 121 says that the noetherian operators are generated by the ones corresponding to 
the corner monomials by taking the closure with respect to the morphisms ([2]). This fact allows 
to come up with a general procedure that constructs the noetherian operators starting with the 
corner monomials and then generates the entire space of noetherian operators. 

Algorithm 3.17. Let I <Z R he a zerodimensional primary ideal of multiplicity fi centered at 
the origin. The following list of instructions construct the noetherian operators associated to I: 

Input; a Grobner Basis G of I and the residual monomials of R/I. 
Output; the space of noetherian operators associated to I. 

• Construct the set C of corner monomials using definition \3.16l 

• For each corner monomial m £ C find the associated noetherian operators by 
rewriting it backwards with respect to Q using definition \3.13l Stop when the backward 
reduction is not possible anymore or when the degree of the polynomial obtained is 
M-1. 

• Collect all the polynomials obtained in the set D. 

• Compute the closure of V by applying the morphism , i = 1 . . . n to all its 
elements. 

• For each element L in the closure ofD calculate D{L). 
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3.4 Extension to modules 



All the results in the previous subsections can be extended in a straightforward fashion to the case 
of zerodimensional primary modules. Rather than giving the details, we use the CoQA version 
of the algorithm for modules to look at a couple of examples. 

Example 3.18. Let A be the matrix 




A 



and let M be the module generated by the rows of A, i.e. M = {xei + 62, xe2 + yei, 1/62). The 
module term ordering we choose is Lex-Pos, meaning that to compare two terms we first look 
at the power product, using Lex, and then we look at the position. The way we just wrote the 
generators of M reflects this choice. It is clear that Jm = {x^ — y, J/^, xy), and, using for example 
CoOcA, we find out that: 

- Ai(M) = 3 

- the Lex-Grobner Basis of M is ^ = {xei + 62, xe2 + yei, ye2, y^ei} 

- a Macaulay basis for M is the set {ei,e2,yei}. 

We begin by writing explicitly the vectorial Taylor expansion of a vector w{x,y) S up to 
degree 2: 

T2w{x,y) = cIqSi + cgoe2 + cl^xei + cfoa;e2 + cj^yei + cl^ye2 + (^^x^ei + (?2^x^^i + Cnxyei + 
Cnxye2 + 42^^^! + 0^22/^62. 

Only few terms survive after we compute the normal form relative to the Grobner Basis Q, 
leading to 

NF,^ (w) = [4o]ei + [4q - cio]e2 + [cIq + cji - cl^jyei. 
We conclude that the noetherian operators associated to M, written in vectorial form, are 

1 



1)^0 = (1,0), D'oo = {-dx,l), D, 



01 



i-dx- 



+ dy, —dx) 



and it is easy to check that they generate a closed subspace since cJa;(-Doo) 



-'^00- 



Z^QQ and 



Example 3.19 (Solution of a system of PDEs). In the introduction we saw that the Funda- 
mental Principle can be used to write an integral representation of the solution of a system of 
linear constant coefficient partial differential equations. We will show how this can be applied, 
now that we know how to compute noetherian operators. Consider the overdetermined PDE 
system given by 

fzz - fz + ft + '^gz = g 

fzt + gt = 

ftt + gzt - gt = 
. ft - gzz + gz + gt = 

where f,g G C°°(M2) and we use indices to denote derivatives. The general solution to ([8]) can 
be written using a generalization of ([1]). We consider the rectangular operator P{D) defined by 



X + y 2a; — 1 \ 



P 



xy 
y 



y 

xy -y 



y J 
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where x and y are the dual variables of z and t respectively. Note that we are choosing a 
particular Fourier transform to write P{D) so that it does not take into account the factor 
—^J—\. The module M associated to the matrix P is not primary, hence we can use Singular to 
get a primary decomposition (using the function modDec form the library mprimdec . lib). M is 
the intersection of the two zerodimensional modules 
Ml = ((x,l),(y,x),(0,y)), J^ = ^^ = {x,y) 
M2 = ((x-l,l),(y,0),(y,x-l)), J2 = = 

of multiplicity, respectively, 3 and 2. We already computed the operators associated to the 
module Mi in the previous example. To compute the operators associated to M2 we need to 
shift the variety to the origin using the change of coordinates {X = x — 1,Y = y). Then, 
using the new variables X and y, we can apply the module version of algorithm 13.81 and find 
the noetherian operators: {(1,0), {dX, —1)}. Going back to the variables x,y we have the set 
{(1,0), (9x, —1)}. Therefore, it is possible to write explicitly the solutions to ([8]) as 



1(0,0) ^ ' 1(0,0) ^ ' 1(0,0) 



/ 1 \ ( dx\ +ty _( A-BZ + \Cz^ + Ct + De' + Eze' 

[Oj. ^ I -1 j, I B-Cz-Ee' 

1(1,0) ^ ^ 1(1,0) ^ 

4 The case of positive dimension 

When dealing with ideals and modules whose dimension is positive, in general one may not 
expect the associated noetherian operators to be constant coefficient linear operators. In fact, 
this is the case for some of the examples from the literature (see [71 HH]). For instance, when 
considering the ideal I = (x^, y^, —xz + y) C C[x, y, z\ one has that a set of noetherian operators 
associated to / is {l,dx + zdy} and it can be proved that there exist no set of noetherian 
operators with constant coefficients associated to / (see |19) . example 4, p. 183). However, an 
interesting property that we notice in this case is that the set of "differential" variables from the 
set of variables appearing in the polynomial coefficients (in this case such sets are respectively 
{x, y} and {z}). This is actually valid whenever we can put the algebraic variety in a particular 
position, through an opportune change of coordinates, called normal position. To do this, one 
can apply the procedure of Noether normalization to the ideal /. This algorithm comes from 
the so-called Noether Normalization Theorem (see [1], p. 116). We now state a version of the 
theorem that we will need for our computations: 

Theorem 4.1 (Noether Normalization Theorem). Let I be a primary ideal of C[zi, . . . , Zn]. 
There exist a non-negative integer d and a (linear) change of coordinates 

if : C[zi,.. .,Zn]^ C[xi,. . .,Xn-d,tl, • • • , ^d] 

such that: 

a) ^{I)nC[ti,...,td] = (0), 

b) C[zi, . . . , Zn]/I is a finitely generated C[ti, ... ,td]-module, 

c) for each i = 1 . . . n — d, (p{I) contains a polynomial of the form 

Qi{ti, . . .,td,Xi) = x1' +pi{ti, . . . ,td)x^*"^ 4 hPe,(ii, • • • ,*d) 

where Cj is the degree of the polynomial Qi . 

The ideal <f{I) is said to be in normal position with respect to the variables xi, . . . ,Xn-d- 
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Remark 4.2. The proof of the Normahzation Theorem can be found for example in pj, in 
the case of prime ideals. However, as shown in [10], the result holds for the general case with 
the exception of condition a) which requires I to be primary. If the ideal / is prime, the 
polynomials Qi in condition c) can be chosen to be irreducible. The proof of the theorem 
provides an algorithm to achieve the normal position. Basically, at each step one constructs the 
polynomial Qi, performing a generic coordinate change such that Qi has a monic leading term 
of the form x^' , and then one eliminates the variable Xi . A procedure to compute the Noether 
normalization of an ideal has also been studied in [13] and it is available in Singular through the 
library algebra, lib (see [9j and its manual). We coded a version of the algorithm for OoCoA as 
well, gj. 

Theorem 14.11 basically states that it is possible to find a new system of coordinates where the 
X variables act as "variables" and the t variables act as "coordinates", and where the integer d 
appearing in 14.11 is nothing but the dimension of the ideal /. Hence, if we make the variables 
t invertible, i.e. if we extend the ideal to the ring C(t)[x] where C(i) is the ring of quotients 
of C[t], we end up with a zerodimensional ideal. Furthermore, since we are interested only in 
primary ideals, we may expect that the extension of the ideal to C(t)[x] is still primary. The 
following proposition assures that such facts hold if / is in normal position. 

Proposition 4.3. Let I = (/i,...,/r) be a primary ideal of dimension d in the polynomial 
ring R = C[xi, . . . , Xn-d, ti, ■ ■ ■ , td], in normal position with respect to xi, . . . , Xn-d- Denote by 
Rd = C{ti, . . . ,td)[xi, . . . ,Xn-d] the ring of polynomials in the x variables with coefficients in 
the field of fractions C{ti, . . . ,td) = Frac(C[ti, . . . , td])- The following facts hold: 

1) the inclusion map ip^^ : I ^ IR^ is infective and IR^ Ci R = I, 

2) the extended ideal IRd is primary, 

3) the extended ideal IR^ is zerodimensional. 

Proof. The fact that the inclusion is injective is trivial. To prove 1), let us consider a polynomial 
/ in i? n IRd- As an element of IR^ it can be written in the form 



where x = {xi, . . . ,Xn-d), t = (ti, ■ ■ ■ ,td), and Oj and hi are just polynomials in the set of 
variables indicated in parenthesis. Let b{t) = ni=i^i(*) consider the product bf. Both b 
and / are polynomials in R and their product is an i2-linear combination of the generators of /, 
so bf E /. Since / is primary it follows that either b"^ €z I for some positive integer m or f G I. 
The first possibility is in contradiction with condition a) of the Noether normalization, hence 
f £ I. This proves that IR^ n i? C /. The opposite inclusion is trivial, so we conclude that 
IRd n R = I. The same type of argument can be used to prove that IRd is primary: consider 
two fractions 

a{x,t) c{x,t) 

such that fg £ IRd- Then (bf) ■ {eg) is a polynomial in / and since / is primary we either have 
6/ G / or e^g"^ G / for some positive integer m. In the first case, using again that / is primary 
and using condition a) of Theorem 14.11 we get that / is in /. In the second case we have that 
5™ is in /. Therefore either / G IRd or 5™ G IRd- Finally, statement 3) follows from the theory 
of the dimension of an ideal, since (ti, . . . , td) is a maximal regular sequence in R/I that reduces 
to just constants when extending the ideal to -R^. □ 
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Before we move on and present an equivalent version of algorithm [3]8] for non zerodimensional 
ideals, there is still one more step. Formerly, when treating the zerodimensional case, we chose 
to start with a Grobner Basis for the ideal /, computed with respect to any term ordering. This 
is no longer possible if we want to extend the procedure to the positive dimensional case. In 
fact, after we perform the normalization, the variables t play the role of "constants" once we 
extend I to Rd = C{t)[x]. The following example illustrates a problem that may occur if we do 
not choose carefully the term ordering on R. 

Example 4.4. Consider the ideal / = (x^ — t,xt — 1) in C[x,t]. A DegLex-Grobner Basis for 
/ (with X > t) is given hy Q = {x^ — t,xt — l,t^ — x}, where the leading term are highlighted 
in bold. When we look at such polynomials in R^, however, we see that the leading terms 
change, in fact the last polynomial should better be written as — x + 1^. Note that in this case 
the extended ideal IR^ happens to be the whole ring R^ since the polynomial — 1 belongs 
to IRd, and such polynomial is a constant in C(t)[x]. It is a necessary and sufficient condition 
for an ideal to be the whole ring that any Grobner Basis with respect to any ordering contains 
a constant polynomial, but if we look at Q we see that there is no such a constant, meaning a 
polynomial only in the variable t. Therefore we conclude that the set Q does not form a Grobner 
Basis for IRd, with respect to the ordering DegLex restricted to the terms in x. If we choose 
instead the term ordering Lex, a Grobner Basis for / is given by ^ = {— x + 1^, — 1}, and in 
this case it contains a polynomial in t, making Q a Grobner Basis for IRd as well. 

As the example shows, we really want the variables x to be the main variables with respect 
to which the Grobner Basis needs to be computed. This can be achieved using Lex, or any 
other elimination ordering with respect to the variables x. Lemma 12.31 then ensures that after 
extending the ideal to Rd, Grobner Bases are preserved. We now have all the ingredients to 
generalize algorithm 13.81 to the case of an ideal of dimension greater than zero. As in section [3l 
we will suppose that a primary decomposition of the ideal has already been calculated. 

Algorithm 4.5 (Noetherian operators for positive dimensional ideals). Let d be a positive inte- 
ger, X = (xi, . . . , Xn-d) CLud t = {ti, . . . ,td) be variables and let a = ax-fTt be a product ordering. 
Let I be a primary ideal in R = C[x,t]. Suppose that I is in normal position with respect to 
X. Moreover, let IRd be the extended ideal in Rd = C{t)[x] and suppose that the characteristic 
variety of IRd inC{t)'^ is the origin. The following procedure computes the noetherian operators 
associated to I: 

Input.' Q = {gi, . . . , gr} a o" -Grobner Basis for I. 
Output; a set of noetherian operators for I. 



• Compute the multiplicity of the ideal, fJ-{I). 

• Write the Taylor expansion at the origin of a polynomial h £ C[x] 
up to the degree — 1 with variable coefficients Ca-' 



a\<fi 




. . .X, 




(9) 



qGN' 



• Let be the leading term of gi and define := t^^ ■ ■ ■ t^^ 

Repeat 
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• Multiply h by t"' and compute its normal form with respect to Q . 

• Rename that as h: 

h := NF, (tU) = ^ (i^(t)xf . . . xfr/ (10) 
Until the number of nonzero dp is exactly ^. 

• For each j3 such that dj^ ^ 0, find polynomials apa{t) such that dj^{t) = 
Ylia'^Pct{t)ca o,nd return the operator 

^P = Y1 — i — ~ 1 "^^i' ■ ■ ■ ^^n-'d = Yl (^l3a{t)D{ai, a^-d, 0, . . . , 0). 

0!i . • ■ ■ Q-n—d- 
a a 

Proof. Let hhe a polynomial of R. We want to characterize the membership of h to I. Since we 
are assuming that / is in normal position, by condition 1) of Proposition 14.31 this is equivalent 
to the membership of h to IRd- Since the latter is a zero dimensional ideal of multiplicity 
h E IRd if and only if the Taylor polynomial of degree fi — 1 of h, with coefficients in C{t), 
reduces to zero when rewriting it using a Grobner Basis for IRd- This follows from the the same 
proof as in algorithm 13.81 By Lemma 12.31 a u^^-Grobner Basis for IRd is given by the same 
elements of the Grobner Basis of /. This means that computing a normal form in / and in IRd 
is equivalent. However, when writing the Taylor expansion ([9]), we need to consider that the 
coefficients Ca also depend on t. In order to be able to perform a one-step reduction, we need 
each term in ([9]) to be at least multiplied by t"'. This does not affect the membership of T^-ih 
as a polynomial in Rd since it is just a multiplication by a constant. Also when considering 
the expression ([9]) in C[x,t], the effect of the multiplication does not change the annihilation of 
NFo- (T^-i/i), since obviously 

NF, {T^-ih) = 4^ NF, {f<T^-ih) = 0. 

The one-step reduction is then iterated enough times in (jlOp until we reach a sufficiently small 
number of nonzero terms (namely fi). By what we have proved so far, it is then clear that at 
the end of the process the polynomial h is exactly the normal form of T^-i/i as a polynomial in 
Rd and hence the annihilation of its coefficients is equivalent to the condition h G IRd- D 

Remark 4.6. The main difference with respect to the algorithm for zerodimensional ideals 
is that, in this case, we do not know if after just one step of reduction we have achieved the 
normal form of the polynomial h{x,t), since the multiplication by could not be enough to 
assure that h has been rewritten to a sum that runs over just the Macaulay basis terms for IRd- 
Multiplying T^_i^ once by is definitely enough for a one-step reduction of each term of the 
Taylor expansion. That is, each term is being rewritten using at most one of the elements of 
the Grobner Basis. However, further reductions might occur if we multiply again by f^. Also, 
note that such an iteration has to terminate because ax is a well ordering. 

Remark 4.7. The reduction step (jlOp for ideals with few generators is not very heavy, but 
performing it multiple times could slow down the procedure by a significant amount. We believe 
that it is possible to find an exponent 71 large enough so that we need to multiply by just 
once, allowing the reduction to bring h all the way down to its final expression. For example, 
choosing 71 = ^-7 seems to work fine at least in the cases we tested, without the need of further 
iteration. 
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When applying algorithm 14.51 to an ideal / in normal position, some redundant factors in t 
could appear as an effect of the iterative multiplication by V at each step. Since such factors are 
constants in Rd, they are actually not needed to characterize the membership of a polynomial 
in Rfi. It is then possible to eliminate these factors from the final expression of the noetherian 
operators. The next example will clarify what we mean. 

Example 4.8. Consider the system of partial differential equations in three variables given by 

fxx — 

fyy = 

fy ~ fxt 

Its solutions are differentiable functions of the form /(x, y, t) = A(t) + B{t)x + B'(t)y, where A 
and B are arbitrary functions of t. We want to derive this last statement using the fundamental 
principle. The primary ideal associated to the system is / = —xt + y) in C[x,y,t] (see 

[19]). If we consider the Lex ordering where x > y > t, a Grobner Basis for / is given by 
y^, —xt + y). Let us compute the associated noetherian operators using algorithm 14.51 
It is immediate to check that / is in normal position with respect to x and y and that, after 
inverting t, the variety associated to IC{t)[x,y] is the origin in C.{tf. The multiplicity of / can 
be computed with CoCoA., and it is = 2. So we just need to write a linear polynomial h with 
variable coefficients and multiply it by i, which is the only term in t appearing in the leading 
terms of the Grobner Basis: 

Tih = t ■ Tih = tcoo + tciox + tcoiy. 

The only rewrite rule that we need to use to reduce h is hence xt ^ y which leads to the final 
expression for the normal form 

NF^ (h) = [teoo] + [cio+teoi]y. 

Since the terms in x and y of the last expression are exactly = 2, we do not need to proceed 
further and then we conclude that the noetherian operators are {t,dx + tdy}. Since the first is 
a multiple of i, we can divide it by t and get the final set {l,dx + tdy}. Now we can write the 
integral formula for the general solution of the system, using ^, rj, r as dual variables: 

f{x,y,t)= J e^(<+s^^+*")d^i(C,r?,r)+ J (^C + rar?)e^«+J"'+*-)(i^2(C, t) = 

Jr Jr Jr Jr Jr 

The last expression gives exactly the general solution as anticipated above. One just has to 
consider arbitrary Radon measures d^i{T) = A{T)dT and dfj,2{T) = B{T)dT where A and B are 
the Fourier transforms of the two arbitrary functions A and B. 
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